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ABSTRACT 

I present an analysis of the deepest X-ray exposure of a radio millisecond pulsar (MSP) to date, an 
X-ray Multi Mirror-Newton European Photon Imaging Camera spectroscopic and timing observation 
of the nearest known MSP, PSR J0437-4715. The timing data clearly reveal a secondary broad X- 
ray pulse offset from the main pulse by ~0.55 in rotational phase. In the context of a model of 
surface thermal emission from the hot polar caps of the neutron star, this can be plausibly explained 
by a magnetic dipole field that is significantly displaced from the stellar center. Such an offset, if 
commonplace in MSPs, has important implications for studies of the pulsar population, high energy 
pulsed emission, and the pulsar contribution to cosmic ray positrons. The continuum emission shows 
evidence for at least three thermal components, with the hottest radiation most likely originating from 
the hot magnetic polar caps and the cooler emission from the bulk of the surface. I present pulse 
phase-resolved X-ray spectroscopy of PSR J0437-4715, which for the first time, properly accounts for 
the system geometry of a radio pulsar. Such an approach is essential for unbiased measurements of the 
temperatures and emission areas of polar cap radiation from pulsars. Detailed modelling of the thermal 
pulses, including relativistic and atmospheric effects, provides a constraint on the redshift-corrected 
neutron star radius of R > 11.1 km (at 3a conf.) for the current radio timing mass measurement of 
1.76 Mq. This limit favors "stiff" equations of state. 

Subject headings: pulsars: general — pulsars: individual (PSR J0437-4715) — stars: neutron - 
X-rays: stars — relativity 
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1. INTRODUCTION 

PSR J0437-4715 was discovered by Uohnston et alj 
(1993) in the Parkes southe rn radio pulsar su rvey. At 
a distance of 156.3 ± 1.3 pc (jDeller et alJl2008f ). it is the 
nearest known rotation-powered "recycled" millisecond 
pulsar (MSPs). It has properties typical of the Galac- 
tic population of MSPs, with a spin period P = 5.76 
ms and spin-down rate (after kinematic corrections) of 
P = dP/dt — 1.0 x 10~ 20 s s -1 , corresponding to a 
dipole magnetic field strength B = 3 x 10 8 G, a char- 
acteristic age r w 4.9 Gyr, and spin-down luminosity 
E — 3.8 x 10 33 ergs s _1 . The pulsar is bound to a 0.2 
Mp) white dw arf companion in a 5.74 day circular orbit 
(|Bailvnlll993ft . 

In soft X-rays, PSR J0437-4715 was the first 
MSP to be firmly detected, in a serendipitous dis- 
cover y in the course of the ROSAT all sky sur- 
vey ([Becker fc Trumpeil 11993ft. Subsequent observa 



tions with ROSAT (|Becker fe Trumped 11999ft . the Ex- 
treme Ultravio l et Ex plorer Deep Survey Instrument 
" T99l . Chandra ACIS-S and HRC-S 
and XMM-Newton EPIC MOS and 



Halp ern et alJ 
Zavlin et alJ Tg002) . 



pn (jZavlinl l2006f ) revealed a single broad, asymmetric 
pulse and a relatively soft spectrum, with at least two 
distinct spectral components. The bulk of emission up 
to ^2 keV is thermal in nature. The small inferred ef- 
fective radii (<2 km) suggest that this radiation is as- 
sociated with the neutron star polar caps. Given the 
large characteristic age of the pulsar, the heat is most 
likely not due to neutron star cooling but is instead 



continuously supplied by an energetic return current 
from the magnetosphere ab ove the polar c aps (see, e.g., 
lHarding fc Muslimovll2002t) . Above ~2 keV. I Zavlin et all 
(2002) have identified an additional faint component that 
is best described by a power-law. The nature of this spec- 
tral tail is less certain as it could be due to a variety of 
plausible mechanisms such as pulsed non-thermal radi- 
ation due to particle acceleration in the pulsar magne- 
tosphere, unpulsed synchrotron radiation due to an un- 
resolved, faint compact pulsar wind nebula or intrabi- 
nary shock, or a spectral tail caused by inverse Compton 
scattering of the soft thermal photons by particles (most 
likely electrons and p ositrons) of small optical depth 
fBogdanov et al. 2006). PSR J0437-4715 has also been 
detec ted in the UV (|Kargaltsev et ani2004tlDurant et al.l 
12012ft and recently identified as a pulsed source in 7- 
rays (>100 MeV) with the Fermi Large Area Telescope 
(jAbdo et al.l[2009ft . 

Based on existing X-ray observations, modeling of 
the pulsed emission from PSR J0437-4715 has revealed 
unique insight into the properties of the stellar surface, 
magnetic field configuration, and neut r on star structure. 
For in stanc e. [Pavlov fc Zavlin! (11997ft . iZavlin fc Pavlovl 
(fl998ft . and iBogdanov et al.l (|2007ft~ have demonstrated 
that a hydrogen atmosphere is likely present at the 
neutron star surfa ce, as expected for a recycled pulsar 
(Alpa r~et al.lll982T ). The morphology of the thermal pul- 
sations seen in the archival XMM-Newton data point to 
an offset of the magnetic axis from the neutron star cen- 
ter (jBogdanov et al. 2003). In addition, the compactness 
of PSR J0437-4715 is found to be R/R s > 1.6 (where 
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Fig. 1.— XMM-Newton MOS1 and MOS2 mosaic image of the field around PSR J0437-4715 in the 0.1-10 keV band. The rectangles 
show the relative positions and orientations of the active EPIC pn chip 4 used in observations 0603460101 {white) and 0112320201 {green), 
with the timing readout direction indicated by the arrow. The dashed lines delineate the EPIC pn detector columns used to extract the 
source counts from the pulsar. Note the significant contribution of photons from the background source RX J0437. 4-4711 in the source 
extraction columns of observation 0112320201. 



R s = 2GM/c 2 is the Schwarzschild radius) at 99.9% 
confidence, which for the best available m ass measure- 
ment (1.76 ± 0.20 iVerbiest et al.ll2008D corresponds 
to R > 8.3 km (iBogdanov et al.ll2007|) . As shown in 
iBogdanov et al.l (|2008|) . in principle, with deeper X-ray 
exposures it may be possible to obtain even tighter lim- 
its on the allowed neutron star equations of state. As 
the nearest and brightest known MSP, PSR J0437-4715 
is the best-suited target for detailed investigation, espe- 
cially given the availability of an independent mass mea- 
surement derived from radio timing observations. This 
combination can potentially provide stringent constraints 
on neutron star structure. 

In this present paper, I report on a new XMM-Newton 
spectroscopic and timing observation of PSR J0437- 
4715, which sheds new light on the X-ray properties of 
MSPs and the physics of neutron stars, in general. The 
paper is organized as follows. In §2, 1 overview the obser- 
vations and data reduction. In §3, I describe the char- 
acteristics of the pulse profile, In §4 and §5, I present 
phase-averaged and phase-resolved spectroscopic analy- 
ses of the pulsar. In §6, I present modelling the pulsed 
emission from PSR J0437-4715. I offer conclusions in §7. 

2. OBSERVATION AND DATA REDUCTION 

PSR J0437-4715 was revisited with XMM-Newton 
between 2009 December 15 and 17 (observation ID 
0603460101) in a continuous 130-kilosecond exposure, 
corresponding to the entire available observing time of 
revolution 1835. The European Photon Imaging Cam- 
era (EPIC) MOS1/2 instruments were configured for full 
imaging mode. The EPIC pn was used in fast timing 
mode, in which only detector chip 4 is active, enabling 30 
/is relative timing precision while sacrificing one imaging 
dimension. For the MOS1/2 and pn, the thin optical fil- 
ter was in place. The dispersed Reflection Grating Spec- 
trometer (RGS) data yielded no interesting spectral or 
timing information, owing to the relatively faint nature 
of the pulsar, and were thus not used in this investiga- 
tion. 



The data reduction,imaging, and timing analyses were 
carried out using SASJ 10.0.0 and FTOOL^I 6.10, while 
the spectral analysis was conducted in XSPEC0 12.6.0q. 
The raw (ODF) MOS and pn datasets was reprocessed 
with the SAS emchain and epchain tools, respectively, 
and filtered for instances of strong background flares. 
Discarding time intervals with severely contaminated 
data results in 104.9, 104.9, and 120.4 ks of usable expo- 
sure time for MOS1, MOS2, and pn, respectively. Subse- 
quently, the recommended processing filters (flag, pulse 
invariant, and pattern) were applied to generate the data 
used in the spectral and pulse profile analyses. 

For the phase-integrated spectroscopy, the counts from 
the pulsar in the MOS 1/2 data were extracted from 60" 
circles cent ered on the source p osition derived from ra- 
dio timing (|Verbiest et al.ll2008f) , which contain ^90% of 
the total energy at ~1.5 keV. The background was es- 
timated by considering three source-free regions in the 
immediate vicinity of the pulsar. For the phase-resolved 
spectroscopy and pulse profile modelling of the EPIC pn 
dataset, the counts from the MSP were obtained from 
RAWX columns 29-43 (inclusive), corresponding to a 
30.75" circle in full imaging mode, which encloses ^88% 
of the total energy for 0.3-2 keV. This narrow extraction 
region was chosen so as to minimize the large background 
level (owing to the one-dimensional imaging mode used 
for the pn detector), which dominates beyond ~30" from 
the source position. To obtain a reliable estimate, the 
background was taken from source-free readout columns 
on both sides of the columns containing the pulsar. For 
the spectroscopic analysis of the EPIC MOS1/2 and pn 
data, the extracted counts were grouped using at least 
30 and 150 counts per bin, respectively. 

The archival XMM-Newton observation of PSR J0437- 

1 The XMM-Newton SAS is developed and maintained by the 
Science Operations Centre at the European Space Astronomy Cen- 
tre and the Survey Science Centre at the University of Leicester. 

2 Available for download at http:/ /heasarc. gsfc.nasa.gov/ftools/ 

3 See http://heasarc.nasa.gov/docs/xanadu/xspec/ 
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Fig. 2. — A comparison between the XMM-Newton EPIC pn X- 
ray pulse profile in the 0.3-2 keV range (middle panel), the tem- 
plate profile at 1.4 GHz (top panel), and the Fermi LAT profile 
(bottom panel) of PSR J0437-4715. Phase zero is defined based on 
the radio ephemeris. The error bar in the middle panel corresponds 
to the 70 fis absolute timing uncertainty of XMM-Newton. The 
vertical dotted lines mark the phase intervals used in the phase- 
resolved spectroscopic analysis discussed in §5. 



4715 (observation ID 0112320201, revolution 519), with 
the same instrument configuration as the new observa- 
tion, was performed with an unfavorable mean telescope 
position angle (137°). As a consequence, the readout 
columns containing the pulsar emission are significantly 
contaminated by t he background active galac tic nucleus 
RX J0437.4-4711 (jHalpern fc Marshalll fl996h . In con- 
trast, the new observation was carried out at a position 
angle that avoids the AGN entirely (192°, see Figure 1) 
thus resulting in significant improvement in the quality 
of the timing data and a more reliable estimate of the 
background in the EPIC pn. In light of the problems 
with the archival data, the two EPIC pn observations 
were not combined and only the new, uncontaminated 
dataset was used. 



3. PULSE PROFILE ANALYSIS 

To obtain the X-ray pulse profile of PSR J0437-4715, 
the arrival times of the photons extracted from the pn 
source region (Fig. 1) were translated to the solar system 
barycenter with the barycen tool in SAS , assuming the 
pulsa r position derived from radio timing ( Verbicst et al. 
120081) and the JPL DE405 solar system ephemeris. The 
corrected arrival times w ere folded at the radio ephemeris 
of PSR J0437-4715 from lVerbiest et all (|2008l) using the 



TEMPO/Q pulsar timing package. For an illustrative 
comparison, Figure 2 shows the radio, X-ray and 7-ray 
pulse profiles of PSR J0437-4715 all folded using the 
same ephemeris, hence showing the absolute phase align- 
ment of the three pulse profiles. 

Previous observations have uncovered that the X- 
ray pulsations from PSR J0437-4715 is characterized 
by a single broad and asymmetric peak that is sig- 
nificantly wider than the radio and 7-ray counterparts 
(jZavlin et all l2002t iZavlinl 12006ft . The substantial im- 
provement in photon statistics of the new data clearly 
reveals a previously unseen "hump" at phases 0.5-0.75. 
For surface polar cap radiation, this feature likely cor- 
responds to the secondary hot spot on the "far" side of 
the neutron star. The implied peak-to-peak separation 
of the two pulses in phase is ~0.55 as measured from the 
more prominent to the fainter pulse, substantially dif- 
ferent from 0.5 as expected from antipodal polar caps. 
This implies hot spots that are not diametrically oppo- 
site, presumably due to a significant offset of the mag- 
netic dipole axis from the center of the star. Note that 
for PSR J0437-4715, the magnitude of the Doppler ef- 
fect that would be induced by the rapid stellar rotation 
is not large enough to cause the apparent asymmetry, 
espe cially when a H atmo sphere is assumed (see Figure 
1 in lBogdanov et al.ll2007f ). 

The pulsed fraction (defined by convention as the 
portion of counts above the pulse minimum) of the 
background-subtracted lightcurves was determined from 
the pulse profile modelling presented in §6 in five energy 
bands: 0.275-0.35, 0.35-0.55, 0.55-0.75, 0.75-1.1, and 
1.1-1.7 keV. This approach yields 32% ± 1%, 35% ± 1%, 
37% ± 1%, 37% ± 1%, and 35% ± 2%, respectively. It has 
been shown that for the range of plausible NS masses 
and radii, isotropic surface thermal emission from a neu- 
tron star surface, such as that due to a blackbody, can- 
not exceed ~33% even from point-like hot spots for the 
plaus i ble range of NS m asses and radii (jPsaltis et al.l 
I2000t iBeloborodovl I2002T) . Nonetheless, the computed 
pulsed fractions are still in agreement with a predom- 
inantly thermal origin of the observed radiation if a 
H atmosphere is present on the neutron star surface. 
This is because the inherently anisotropic pattern of the 
emergent radiation from such an atmosphere can pro- 
duce much larger rotation-induced modula tions (see, e.g., 
IZavlin et al.|[l99l IZavhn fc Pavlovl [l998). The slightly 
lower pulsed fraction at energies below ^0.35 keV can be 
attributed to the contribution of emission from a larger 
portion of the stellar surf ace, as suggested by the UV 
emission from the pulsar (|Durant et al.l 120121) and the 
spectroscopic analysis in §4, which would be less modu- 
lated due to its greater physical extent. The pulse profiles 
for the other four energy bands are consistent with hav- 
ing identical pulsed fractions. Although the number of 
source photons is quite sparse above 2 keV, faint mod- 
ulations similar to those at lower energies are seen in 
the 2-6 keV band. Closer inspection of the model spec- 
trum shown in Figures 3, however, reveals a significant 
contribution from the hot thermal component, implying 
that the modulation may still be due to thermal emis- 
sion. Unfortunately, due to the high background of the 



4 Available for downl oad 

http : / / www. atnf . csiro.au / research / pulsar / tempo / 
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Fig. 3.— The XMM-Newton X-ray EPIC MOS1/2 spectra of 
PSR J0437-4715 fitted with a two-temperature H atmosphere 
model plus a power-law and a cool blackbody (solid line). The 
middle six panels show the fit residuals for the various models listed 
in Table 1. The bottom panel shows the best fit two- temperature 
H atmosphere model plus a power-law and a cool blackbody. The 
dashed and dotted lines show the individual absorbed and unab- 
sorbed model components, respectively, while the solid line shows 
the total spectrum. See text and Table 1 for best fit parameters 
and acronym definitions. 

EPIC pn data, combined with the significant decline in 
the telescope effective area at higher energies, no use- 
ful pulse profile information is available above ^3 keV, 
where the contribution from the thermal components is 
negligible. Thus, only an uninteresting limit of <90% 
pulsed fraction for the background-subtracted emission 
can be placed, precluding any insight into the origin of 
the power-law tail. 

4. PHASE-AVERAGED SPECTROSCOPY 

iZavlin et all (|2002l ) and lZavlinl (|2006l ) have found that 
the X-ray spectrum of PSR J0437-4715 cannot be ade- 
quately represented by a single or even most two compo- 



nent emission models. Based on this, in the spectral fits I 
apply two plausible multi-component models: (i) a multi- 
temperature thermal plus single non-thermal model, and 
(ii) a multi-temperature Comptonized thermal model. 
For the thermal components, I consider both a sim- 
ple bl ackbody and a non-magnetic H atmosphere spec- 
trum (jRomanill 19871: IZavlin et al.lll996tlMcClintock et al.l 
120041 ). which is appropriate for PSR J0437-4715 given its 
low magnetic field (B « 3 x 10 8 G). Even though polar 
cap heating in MSPs is believed to be due to a return cur- 
rent of relativistic efrom from the magnetosphere above 
the polar caps, the H atmosphere model assumption of a 
heat source beneath the atmosphere is still fully valid. 
This is because the penetration depth of the imping- 
ing relativistic particles from the magnetosphere is much 
greater (~2-3 orders of magnitude) than the characteris- 
tic depth of the atmosphere, meaning t hat their energy 
is dep osited below the atmosphere (see iBogdanov et al.l 
2007, for details). 

Blackbody fits to MSP spectra tend to produce in- 
ferred emission areas that are substantially smaller (a 
factor of ~10) than the classical polar cap area, R pc = 
(2nR/cP) 1 ^ 2 R. In contr ast, a H atmosphere yields 
values comparable to Rr, r (iBecker fc; Asdienbachl [20021 : 
IZavlin! 12009 : IBogdanov fc Grindlavl 12009ft . It has been 
suggested that the discrepancy for the blackbody model 
could be du e to heating of only a small portion of the 
polar caps (|Zhang fe Chengj 120031 ). More plausably, it 
suggests that an atmospheric layer is present on the stel- 
lar surface, as expected from evolutionary arguments. In 
particular, the peak intensity of a H atmosphere contin- 
uum occurs at higher energies than a blackbody for the 
same effective temperature. As a result, fitting a black- 
body to H atmosphere emission would produce a higher 
inferred temperature and, as a result, a significantly 
smaller effective area. Furthermore, blackbody radiation 
has been shown to provide an inadequ ate description of 
the thermal pulsations from this MSP (jZavlin fe Pavlov! 
1 1998t IBogdanov et aLlfeOOTT ). Nevertheless, for complete- 
ness and to enable convenient comparison with previ- 
ously published X-ray observations of MSPs, I have con- 
sidered a blackbody in the spectroscopic analysis as well. 

It is important to emphasize that even for phase- 
integrated spectroscopy it is essential to properly ac- 
count for the rotation of the MSP and the geometric 
configuration of the hot spots, in order to obtain reliable 
measurements of t h e temperatures and emission radii 
dPsaltis et all 120001: IZavlin et all 120021: IBogdanov et~aT1 
l2007f ). For this reason, I consider the model of 
a rotating neutron star with t wo hot spots , base d 
on the prescription describe d in iBeloborodovi ((2 002) . 
iPoutanen fe Gierlihskil (|2003l ). and lViironen fe Poutanenl 
(2004) but with the additio n of a realistic H atmosphere 
(see IBogdanov etldl 120071 for further details). This 
model incorporates the angles between the spin and mag- 
netic axes (a) and the spin axis and the line of sight (£) 
as well as the mass and radius of the pulsar as parame- 
ters in the spectral fits. Thus, the derived effective areas 
are deprojected and redshift-corrected. Moreover, the in- 
ferred temperatures are as measured at the NS surface by 
correcting them for gravitational redshift. For each polar 
cap I consider two concentric emission regions - a small 
high-temperature hot spot surroundend by a cooler an- 
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TABLE 1 

XMM-Newton EPIC MOS1/2 best-fit spectral models and unabsorbed fluxes for PSR J0437-4715 



ModeP 


7V H 
(10 19 cm- 2 ) 


(10 6 K) 


^eff,l b 

(km) 


(10 6 K) 


Rcti,2 b 
(km) 






r 


(0.1-10 keV) 


Xl/dof 


BB(x2)+PL 
Hatm(x2)+PL 


5 7+ 1 - 7 
°-'-1.5 

<0.09 
7 

4.5±3-» 
7 


o n2 +0.14 
J - Uz -0.12 
9 7R+0-03 
z -' D -0.05 

'-'• u ' : '-0.12 

2 14+ - 16 
z - 14 -0.11 
r, nn+0.05 
z - zz -0.04 

2.03lg;l? 


05+ 03 
u - uo -0.02 

u - u '-0.01 

u - u4 -0.02 

n i5+ - 10 

u - lo -0.06 

n io+o.04 

U - iS -0.05 
i 7 +0.05 


1 70+ - 18 

u - yy -0.04 
i 77 +0.15 
'-0.13 
n 75+O.O8 
u - '°-0.08 

u - 00 -o.oi 

n 7^+0.08 


°-n±° ° l 

0.46±g;il 

in+ 05 

u - lu -0.04 

1 15+ 110 

n 65+°- 19 

U.DO_ 25 






9 70 + 0. 04 
z - '8-0.07 

1.56 

9 Qi+0.03 
z - ol -0.03 

2 7n+°- 10 

z - ' u -0.24 

1.56 

9 on+0.04 
z - 8U -0.03 


1.80(1.13) 
1.34(0.20) 
1.88(1.25) 
1.66(0.90) 
1.41(0.15) 
1.85(1.20) 


1.13/282 
1.48/283 
1.13/283 
1.15/282 
1.14/283 
1.15/283 















T cB,3 
(10 6 K) 


^c£F,3 
(km) 








BB(x3)+PL 

Hatm( X 2) 
+BB+PL 


7 

16.7 + |-| 

is.9±r 6 

7 


d - XD -0.10 
o 91 +0.09 
■ 5 " al -0.09 
o 9n +0.08 
o - zu -0.02 

2 56+ - 06 
Z - OD -0.08 
9 C9+0.09 
z - oz -0.09 

2 56+ - 06 
Z - OD -0.09 


05 +U - 02 
u - uo -0.01 

04+ 02 
U - U4 -0.01 

04+ - 02 
U - U4 -0.01 
17 +0.03 
u - 1 '-0.04 

16+ - 03 
u - 1D -0.03 


1 51 +0.08 

1 54+0-06 

1 c-7+0.09 
-0.06 
1 nf) +0.03 
i ' uu -0.12 
n qq+0-06 
u - yy -0 .22 
1 n4 +0.12 


21 +U - ii 
u - zl -o.io 

21+0- 12 
u - zl -0.08 

iq+0- 07 
u - ly -0.08 
1 o+l.O 
L -°-0.6 
1 O+0.5 
1 -°-0.7 
1 i+l.l 
1 - i -0.6 


44+0-00 
u - 4 +-0.04 

46+ 02 
u - 4o -0.04 

5S+ 04 
U-OJ-0.02 

42+ 06 
u - 4z -0.05 

u - 41 -0.04 
n 40+O.O2 
U - 4S -0.06 


4 6 +4 - 8 

4 Q+3.5 
4 - u -3.0 

1.9t ; 9 

5 0+6.1 

9 i+2.2 
z - i -1.7 


1 70+U.35 
1- '8-0.17 

1.56 

1 70+0.35 
J --'" 5 -0.12 
1 19+O.43 
1 - ±z -0.31 

1.56 

, 15+0.43 

1 - lo -0.32 


2.65(0.19) 
2.46(0.16) 
1.69(0.18) 
3.54(0.11) 
3.58(0.13) 
1.71(0.11) 


1.08/280 
1.07/281 
1.08/281 
1.10/280 
1.10/281 
1.10/281 














(10 6 K) 


^?cff,3 
(km) 


T 






CBB(x2) 
CBB(x2)+BB 


<0.1 
I9.0± 7 ;g 
7 


9 fi o+0.U3 
z - tK >-0.05 
o nl +0.09 

3 03+ 11 
• J - u " 5 -0.08 


omtm 

05+ 02 
u - uo -0.02 

05+ - 02 
u - uo -0.02 


Q3+ u uy 
u - yc> -0.03 

1 45+O.O8 

1 4 q+0-l0 


u - oo -0.14 

25+ 13 
u - zo -0.12 
n 99+O.O9 

u - zz -o.io 


0.43+g-g 6 

52+ - 04 
u - oz -0.03 


5.Q+1I 

1 nr+1.05 
i ' yo -1.07 


Q 115+O.UO6 
u - llo -0.005 

102+ 007 
u - luz -0.008 

n inq+o- 007 

U.1UM_q Q07 


1.33 
2.74 
1.68 


1.31/283 
1.07/281 
1.08/281 



a PL is a powerlaw, BB a blackbody, CBB a Comptonized blackbody, and Hatm a non-magnetic H atmosphere model. All uncertainties and 
limits quoted are at a 1<t confidence level. 

b R e ff calculated assuming a distance of 156.3 pc. For the H atmosphere model, the numbers quoted represent deprojected and redshift-corrected 
effective radii of one hot spot on a neutron star with M = 1.76 Mq, R = 13.5 km, a = 36°, £ = 42.4°, Aa = —25°, and A<f> = —20° (see text for 
definition of these parameters). The third, lowest temperature thermal component in all instances is a blackbody. 

c Unabsorbed X-ray flux (0.1-10 keV) in units of 10 — 13 ergs cm -2 s _1 . The values in parentheses represent the flux contribution of the powerlaw 
component, where applicable. 



nular region, as im plied by previous studies (|Zavlin et al.l 
l2002t IZavlinl 120 06). In all instances, to calculate effec- 
tive radii, I consid er the parallax distance of 156.3 pc 
(jDeller et al.l 1200 8). For the H atmosphere model, I as- 
sume a NS with M = 1.76 M Q , R = 13.5 km. I also 
fix a = 36° and £ = 42° and account for the appear- 
ent offset of the secondary hot spot with the parame- 
ters Aa = —25° and A<j) = —20°, as deduced from the 
pulse profile fits (see §6). The various spectral models, 
their best fit parameters and derived unabsorbed fluxes 
are summarized in Table 1. All uncertainties quoted are 
giv en at a la confiden ce level. Following the analysis 
by iDurant et al.l (|2012f) , the fits are performed with all 
parameters free, as well as with a fixed -/Vh = 7 x 10 19 
cm" 2 or a fixed T = 1.56. The latter was found by 
IDurant et al. (2012) to provide the best fit to the broad- 
band non-thermal spectrum of this pulsar. 

Regardless of the spectral model employed, fitting the 
phase-averaged MOS and pn datasets jointly results in 
statistically poor fits, with typical \v ^ 1-3- In con- 
trast, the fits to either the MOS 1/2 or pn data give 
Xv ^ 1-1 for the same assumed models. This points 
to a possible discrepancy in the cross-calibration of the 
detectors especially for the fast timing mode of the EPIC 
pn instrument. Similar issues have been f ound in other 
studies using the pn fast timing mode (jZavlinl 120061 : 
lArchibald etall 12010ft . As a result, only the MOS 1/2 
data were considered for the phase-averaged fits as they 
span a much wider energy range than the pn timing data. 
The MOS1/2 data exhibit significant narrow residuals at 
~1.7 keV . This feature can be attributed to the Si-K 



edgeQ. I gnoring this interval results in significant im- 
provement in the quality of the fits. 

4.1. Multi- Temperature Thermal Plus Powerlaw 
Spectrum 

Based on previous Ch a ndra an d XMM - Newt on ob- 
servations, iZaylin et al.l (|2002D . IZavlinl (|2006l ). and 
IDurant et al.l ( 2012|) have reported that the phase- 
integrated X-ray radiation from PSR J0437-4715 is best 
represented by a two-temperature thermal plus a single 
faint powerlaw components. When applied to the new 
EPIC MOS1/2 data, the same model is only marginally 
acceptable and shows significant broad residuals at en- 
ergies below ~0.8 keV (middle panels in Figure 3), 
similar to those o bserved in ROSAT PSPC data (see 
iZavlin et al.l [2002). This suggest the presence of an 
additional cool thermal component. The atmosphere 
model used in this study is only calculated down to 
T c s — 3 x 10 5 K since below this temperature par- 
tial ionization becomes non-neglig ible ()Zavlin et alfl 996; 
iMcClintock "elaTl [20041 . As a consequence, the third, 
cool thermal component was represented by a blackbody. 
This addition results in a significant improvement in the 
fi t quality. 

IDurant et al.l ()2012f ) have found that the spectral shape 
in the UV is consistent with thermal emission from the 
bulk of the NS surface with an effective temperature 
(1.5 - 3.5) x 10 5 . The MOS 1/2 fits suggest values in 
agreement with the upper bound of this range, although 

5 See http : //xmm2 . esac . esa. int /docs/document s/CAL-TN-00 18 .pdf 
for details. 
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TABLE 2 

EPIC PN RESULTS OF PHASE-RESOLVED SPECTROSCOPY FOR PSR J0437-4715. 



Ah Teff.l Rcff,i b T eSj2 i? eft , 2 b T F x c x'/dof 

Model a (cm^ 2 ) (10 6 K) (km) (10 6 K) (km) (0.5-3 keV) 

BB(x2)+PL <12.5 2.90lo:l? 0.07±g;gi} 1 - 2i -°om - 37 ±o.i8 -°- 14 ti'.77 9.4/6.0/4.9/6.4 1.02/480 

<7.9 2.82±g;^| 0.07t°o| 1.22±g;g| O^t ," 1.56 9.3/6.0/4.9/6.4 1.02/481 

7 2.90l°;gg 0.07±g;g| E23±g^| 0.39±g;j;| — 0.11±^;fl 9.6/6.1/5.0/6.6 1.02/481 

Hatm(x2)+PL <4~3 2.22+^;^ aifff^ oT55^^ 6 -°±2.4 2.35toio 9.3/6.1/4.9/6.2 1.06/486 

<1.5 2.161o;q7 0.18t°'og 0.58lg;°5 5.1±3g 1-56 9.5/6.1/4.8/6.2 1.07/487 

7 2.211°;°* 0.16lo;°g 0-57lo.ol 5 - 5 i3.5 2 - 39 tcu8 9.6/6.3/5.0/6.3 1.06/487 
a Model and parameter definitions are the same as in Table 1. 

b For the blackbody fits, the quoted effective radius is that measured at pulse maximum (corresponding to phase 
interval I in Figure 2). For the H atmosphere, the value is the true, deprojected and redshift-corrcctcd effective radius. 
c Unabsorbed X-ray fluxes (0.5-3 keV) in units of 10 — 13 ergs cm -2 s — 1 for the four phase bins. 



the generally higher values may also indicate a temper- 
ature gradient across the stellar surface, possibly as a 
function of distance from the heated polar caps. 

4.2. Comptonized Thermal Spectrum 

As suggested in lBogdanov et all J2006), the power-law 
tail observed in the spectrum of PSR J0437-4715 may 
arise due to inverse Compton scattering (ICS) of the 
thermal X-ray radiation by relativistic e^ of small op- 
tical depth. It is plausible to assume that such particles 
may be located in the magnetosphere above the polar 
caps. For the spectroscopic analysis of PSR J0437-4715, 
I use the compbb Compt onized blackbody model in XSPEC 
(|Nishimura et al.l [l986). As additional parameters, this 
model includes a scattering particle temperature kT e and 
optical depth r. For the two thermal components the 
t parameter was linked as, presumably, both undergo 
scattering by the same population. The fit of a two- 
component Comptonized blackbody also exhibits signifi- 
cant residuals at lower energies, which can be accounted 
for by the additon of a cool blackbody component. 

As expected, kT e and r are strongly correlated, result- 
ing in acceptable fits over a wide range of the two param- 
eters. Given that we expect kT e to be much larger than 
kT, the particular choice of kT e and the details of the en- 
ergy distribution of the scattering e^- have little impact 
on the outcome of the X-ray spectral fits. Note that the 
inferred optical depth requires a particle density in the 
magnetosphere (from the surface to the light cylinder 275 
km above the surface) above the polar caps to be ~10 4-6 
times the Goldreich-Julian density. It is interesting to 
note that this value is comparable to that required for 
PSR J0737-3039A to explai n the observed radio eclipses 
in th e double pulsar system (Lyutikov 20QH lArons et al.l 
[20051) . 

5. PHASE-RESOLVED SPECTROSCOPY 

The sizable collection of source photons in the pn 
dataset permits an analysis of the pulsar spectrum as a 
function of pulse phase. To this end, I have divided the 
pn photons in the 0.5-3 keV range into four phase inter- 
vals, with the first interval centered on the pulse peak 
(see Figure 2). Due to the overwhealming background, 
the pn data above ~3 keV was excluded from the anal- 
ysis. Moreover, as recommended by the XMM-Newton 
Science Support Center, photons below 0.5 keV were ig- 
nored in order to remove electronic noise artifacts seen 




0.6 0.7 0.8 0.9 1 

Photon energy (keV) 



Fig. 4. — The XMM-Newton X-ray EPIC pn phase-resolved spec- 
tra of PSR J0437-4715 fitted with a H atmosphere model {solid 
lines). The phase intervals used are shown in Figure 3. The bot- 
tom panel shows the fit residuals. The Roman numerals correspond 
to the phase intervals shown in Figure 2. 



in fast timing mode. Above this energy, the importance 
of the cool thermal component seen in the MOS spectra 
is neglegible (as evident from the bottom panel of Fig- 
ure 3). Therefore, only two thermal components plus the 
non-thermal component are used. This choice of energy 
cut also has the consequence of making the fits nearly in- 
sensitive to ./Vh- As evident from Table 2, the power- law 
component is also poorly constrained since the 0.5 — 3 
keV band is strongly dominated by thermal emission. 

The spectral fits were performed on all four spectra si- 
multaneously. For the blackbody fits, the spectra were 
fit with the same temperature, while allowing for the flux 
normalization to vary independently for each spectrum. 
In the case of the H atmosphere, all parameters, includ- 
ing the flux normalization, were fitted jointly among the 
four spectra since the pulse phase dependence of the ge- 
ometric projection of the emission area is taken into ac- 
count internally by the model for a given set of M/R, 
a, and C- In other words, the flux normalization in this 
case represents the true, deprojected area (as measured 
at the stellar surface), which is the same at all phases. 
This was accomplished by taking advantage of a func- 
tionality in XSPEC that permits simultaneous fitting of 
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phase-resolved spectra^. In addition to the anisotropic 
emission pattern, an important distinction of the H at- 
mosphere model from a blackbody is the angular depen- 
dence of the emergent spectrum. In particular, the con- 
tinuum emission becomes softer for larger angles with 
respect to the surface normal as a consequence of the 
temperature gradient within the atmosphere. Observa- 
tionally, this property should appear as a shift in the 
pea k of the spectrum to lo wer energies (see, e.g., Figure 
2 in IBogdano v et al.l 120071 ). Therefore, it is essential to 
properly account for the geometry to ensure an accurate 
measurement of the effective temperature of the radia- 
tion. Furthermore, since the emission area is observed in 
projection it is important to correct for this effect. 

The results of the phase-resolved spectroscopic analy- 
sis are summarized in Table 2. I t is important to note 
that (as shown in prev ious studies IZavlin fc Pavlovlll998t 
iBogdanov et al.l 120071 ) , although the blackbody model 
can reproduce the spectral shapes of this pulsar's phase- 
resolved spectra for the same temperature, the implied 
variations in the effective area cannot be explained by a 
rotating hot spot model. This lends further support for 
anisotropic thermal emission, which, in turn, implies the 
presence of an atmosphere on the surface of PSR J0437- 
4715. 

6. CONSTRAINTS ON THE NEUTRON STAR 
COMPACTNESS 

The observed pulsations of the surface thermal emis- 
sion from PSR J0437-4715 contain valuable information 
regarding the NS compactness, which can be extracted 
via a constraint on the stellar mass-ra dius relation. To 
this e nd, I consider the model from IBogdano v et al.1 
(2007), which assumes a rotating compact NS with 
two identical hot spots, each corresponding to one 
of the magnetic polar caps. As in the spectroscopic 
analysis, each polar cap is assumed to consist of a 
small high-temperature hot spot surroundend by a 
cooler annular region. The model considers a non- 
rotating Schwarzschild metric and incorporates the 
relativistic Doppler effect and pro pagation time delay s 
using the prescription de scribe d i n iBeloborodovl ( 2002). 
iPoutanen fc Gierlihskil (120031). iViironen fc Poutanenl 
(|2004D . and IPoutanen fc Beloborodovl (|2006l ). This 
appr oach is remarkably accurate as long a s P > 3 
ms (jCadeau et al.l I2007L iMorsink et al.1 120071 ). making 
it applicable to PSR J0437-4715. The NS surface is 
assumed to be cover ed by a non-magnetic, optically- 
thick H atmosphere dRomanil 119871; IZavlin et all 119961: 
IMcClintock et al.l \200W . To allow direct comparison 
with the data, the model was convolved with the EPIC 
pn plus thin filter effective area and corrected for the 
encircled energy fraction and the background level. 

The X-ray pulse profile was fitted by considering ten 
parameters: the two temperatures and effective radii of 
each hot spot (Ti, T2, Ri, and R2), the angles a and 
£, the NS radius R, the displacement (in degrees) of the 
secondary hot spot from the expected antipodal position 
(Aa and A</)), and the phase <fi of the main pulse. High- 
precision radio timing observations have measured the 
orbital inclination of the J0437-4715 binary system to 

6 See http://heasarc.nasa.gov/xanadu/xspec/maiiual/ for fur- 
ther details. 
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Fig. 5. — (Top panel) The XMM-Newton EPIC pn pulse profile 
in the 0.5—1.8 keV range with the best fit model (solid line). The 
individual contributions from each polar cap are shown with the 
dashed lines. The dotted line shows the background level. (Bottom 
four panels) XMM-Newton EPIC pn pulse profiles of PSR J0437- 
4715 in the 0.5-0.65, 0.65-0.85, 0.85-1.2, and 1.2-1.8 keV bands 
(from top to bottom, respectively) fitted with a model of a rotating 
neutron star with two-tcmpcraturc H atmosphere polar caps. See 
text for best fit parameters. 



be i = 137.6° (jVerbiest et al.ll2008l ). It is very likely that 
the rotation axis of the MSP has aligned with the orbital 
angular momentum vector as the NS was spun-up during 
the low-mass X-ray binary stage of its evolution. This 
scenario implies £~ 137.6° or, equivalently, C ~ 42.4°. 
Therefore, I fix the value of £ to the latter value. The 
H absorption column along the line of sight was set to 
A^h = 7x 10 19 cm~ 2 , w hile the distance was fixed at 156.3 
pc (|Deller et al.ll2008() . For all intents and purposes the 
0.8% uncertainty in the parallax distance is completely 
negligible. Varying TVh over the range (1 — 10) x 10 19 
cm~ 2 results in virtually no change in the fit since its 
effect is unimportant above 0.5 keV. In the analysis, I 
assumed a fixed mass and allowed R to vary. The proce- 
dure was repeated for M = 0.7, 1.0, 1.3, 1.76 and 2.2 M Q . 
The fit was conducted simultaneously for four photon en- 
ergy ranges (0.5 - 0.65 keV, 0.65 - 0.85 keV, 0.85 - 1.2 
keV, and 1.2 — 1.8 keV) in order to weaken any covari- 
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Fig. 6. — The mass-radius plane for neutron stars showing the 1, 
2, and 3<r confidence contours (yellow, green, and blue hatched re- 
gions, respectively) for PSR J0437-4715. The solid lines are repre- 
sentative theoretical model tracks (from Lattimer & Prakash 2004). 
The horizontal lines show the pulsar mass measurement from ra- 
dio timing (dashe d line) and the a ssociated la uncertainties (dotted 
lines) from lVerbiest et all 12003) . 



ance between the parameters that define the emission 
spectrum (temperatures and effcctve radii) from the ge- 
ometric and relativistic parameters, which do not depend 
on energy. The energy ranges were chosen in such a way 
as to yield comparable numbers of photons for all bands 
and to reduce the contribution of background photons. 
To include the non-thermal component seen in the total 
spectrum, I absorbed the expected number of photons 
per bin into the uncertainties, assuming an unpulsed con- 
tribution and a r — 1.56 powerlaw spectrum. The best 
fit values for each parameter were determined by search- 
ing the x 2 space. To ensure that the these values are the 
true absolute minima of the parameter space, the proce- 
dure was repeated for a wide range of initial values for 
each parameter. Due to the complicated topology of the 
parameter space, the confidence intervals of the best fit 
parameters were determined by way of Monte Carlo sim- 
ulations, using 5 x 10 3 realizations for each of the four 
assumed NS masses. 

The best fits to the X-ray pulsations for the H atmo- 
sphere hot spot model are shown in Figure 5. The plot- 
ted best fit yields xt — 0-9 f° r 65 degrees of freedom. As 
seen previously, unlike a blackbody model, H atmosphere 
emission can account for the morphology of the X-ray 
pulsations and large amplitu de remarkably well. In com- 
parison with the modeling in Bog danov et all (|2007|) . the 
best fit shown in Figure 5 confirms the upward inflection 
in the lightcurve around phases 0.5 — 0.75, which appears 
to arise due to the secondary hot spot. Figure 6 shows 
the M — R plane with re presentative models for diff erent 
varieties of NS EOS (see lLattimer fc Prakashll200lD and 
the limits derived for PSR J0437-4715, with the hatched 
areas representing the 1, 2, and 3er confidence regions. 
The NS radius is constrained to be > 11.1 km at 3a 



confidence, assuming M = 1.76 M Q for all combinations 
of the other parameters. If taken in combination with 
the best available m ass measurement from radio timing 
(jVerbiest et alJ[2~008| ) , the limits are inconsistent with all 
but the stiffest equations of state. It is interesting to 
note that the lower bounds on the confidence intervals 
are much tighter. This occurs because for more compact 
stars it is not possible to produce the large observed pulse 
amplitude, owing to the stronger gravitational bending 
of light effect, which acts to surpress rotation-induced 
modulations. In contrast, the upper bounds allow for 
very large stars since large-amplitude pulsations can be 
produced more easily for a wide range of viewing geome- 
tries. 

The fits to the X-ray pulsations also yield constraints 
on the magnetic field configuration of PSR J0437-4715. 
In particular, assuming M — 1.76 Mg, the magnetic 
inclination is constrainted to be a = 36^3 degrees (at 
la confidence), suggesting a significant misalignment be- 
tween the magnetic and rotation axefl The latitudi- 
nal offet of the secondary polar cap is constrained to be 
Aa = —25^4 degrees. As defined, a negative Aa cor- 
responds to a northward shift from the antipodal posi- 
tion. The displacement in the longitudinal direction is 
found to be A(f> = — 20^ 2 degrees, indicating that the 
secondary spot trails the antipodal position with respect 
to the sense of rotation, as evident from the phase lag of 
the seco ndary pulse in F i gure 5 . Based on Equations 2 
and 3 in iBogdanov et a l. (2008) and assuming the nom- 
inal best fit NS radius of 13.5 km, the corresponding net 
offset of the secondary hot spot from the antipodal po- 
sition across the NS surface is ~6 km, while the implied 
displacement of the magnetic dipole axis from the stellar 
center is Ax ~ 3 km. 

7. CONCLUSION 

In this paper, I have presented an analysis of the deep- 
est X-ray observation of a radio MSP to date, namely, an 
XMM-Newton exposure of PSR J0437-4715. The data 
permit a significantly better characterization of the spec- 
tral continuum of this pulsar in the 0.1-10 keV band 
compared to previous studies. Although the observed 
radiation from J0437-4715 can be described well by mul- 
tiple models, in all cases the vast majority of soft pho- 
tons (>90%) appear to be thermal in nature. This ther- 
mal radiation exhibits at least three components, with 
the hottest two having total effective areas consistent 
with the expected polar cap size. The coolest com- 
ponent, on the other hand, appears to cover a signifi- 
cant portion of the stellar surfac e, consistent with find- 
ings based on UV observations (|Kargaltsev et all 120041 : 
iDurant et al.1 12012). Since passive cooling of any pri- 
mordial heat stored in the NS cannot account for the 
high global NS temperature (>10 5 K), various heat- 
ing mechanisms have been envoked to explain the ex- 
cess, including: re-radiation of the energy deposited 
in the inner crust by polar cap heating, rotochemi- 
cal heating (|Fernandez &; Reiseneggerl [200 5;) . "gravito- 
chemical" heating caused by a time-varying gravitational 
constant (|Jofre et al.1 [2006) , or capture of exotic par- 

7 By convention, a is reckoned from the spin pole towards the 
equator. 
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tides in the NS interi or (lHannestad et all 12 002). See 
iKargaltsev et all (|2004f ) and lDurant et all (|2012l )" for fur- 
ther discussions of the various possible heating mecha- 
nisms that may operate in PSR J0437-4715. 

An important goal of future observations will be to es- 
tablish the true nature of the hard photons (> 3 keV) 
from PSR J0437-4715. This can be best accomplished 
through detailed phase-resolved spectroscopy. Better 
characterization of the optical and 7-ray emission from 
this pulsar may offer additional clues into the nature of 
this radiation. In addition, with the recent launch of 
NuSTABQ a hard X -ray focusing telescope covering the 
6-79 keV range, further insight could be gained into the 
connection between the harder X-ray and 7-ray emission 
from this pulsar. 

For the first time, I conduct phase-resolved X-ray 
spectroscopy of a radio pulsar that correctly accounts 
for the system geometry and the radiation pattern of 
the atmosphere by jointly fitting the spectra from all 
phase intervals. With greatly increased photon statis- 
tics, this potentially powerful approach can yield strin- 
gent constraints on the NS parameters by utilizing a 
two-dimensional approach to spectroscopy in the photon 
energy-pulse phase plane. 

The model pulse profiles in Figure 5 imply a signifi- 
cant offset of the secondary hot spot from the antipodal 
position, which in turn, suggest deviation from the stan- 
dard cente red magnetic dipole configur ation. As demon- 
strated bv lHarding" &: Muslimov (20 ll|), even a slight az- 
imuthally asymmetric displacement of the pulsar mag- 
netic field can substantially amplify the accelerating elec- 
tric field on one side of the polar cap, which when com- 
bined with a smaller field line radius of curvature, leads 
to larger pair multiplicity. As a consequence, the death 
line for producing by curvature radiation shifts down- 
ward in the period-period derivative diagram, resulting 
in a larger number of pulsars having high pair multiplic- 
ities. If common among MSPs, and pulsars in general, 
these magnetic field distortions could have profound im- 
plications for studies of the pulsar population, modelling 
high energy pulsed radiation, and constraining the pulsar 
contribution to cosmic ray positrons. 

The analysis presented herein confirms that the X-ray 
pulsed fraction od PSR J0437-4715 is consistent with 



the existence of an atmospheric layer on the NS sur- 
face and cannot be reproduced by a blackbody. Mod- 
elling of the thermal pulsations gives a limit on the 
NS radius of > 11.1 km (3cr conf.) assuming a mass 
of 1.76 Mq. This lower bound is much stricter com- 
pared to the limit obtained from the shallower archival 
observations. This value is in agreement with several 
NS radius mea surements using quies cent low-mass X-ray 
binaries (e.g., iRutledge et alj 120011: iHeinke et all 120061: 
iWebb k Barretl 120071 ) and X-ray bursts (see, e.g., lOzell 
120061 : lOzel et al.ll2010l ). However, it is important to em- 
phasize that these methods have resulted in a wide range 
of often contradictory constraints on the NS radius (see, 
e.g.. lSuleimanov et al.ll201 fllQzel et alll2012l for the case 
of X-ray bursts) . This highlights the clear need to investi- 
gate the systematic uncertainties and reliability of these 
model-dependent approaches to constraining NS struc- 
ture, including the method employed in this present pa- 
per. 

This issue n o twith standing, as demonstrated in 
iBogdanov et al.l (|2008l ). deeper X-ray observations of 
PSR J0437-4715 should result in even stricter constraints 
on M/R, especially when combined with further refine- 
ment in the NS mass measurement from radio timing 
and joint mod elling of the pulsed X-ray and 7-ray emis- 
sion (see, e.g.. lVenter et all l2009) to better constrain the 
system geometry. Such efforts may produce definitive 
measurements of the M — R relation, and by extension, 
the state of cold, dense matter. 
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